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ABSTRACT 

We exhibit static solutions of multi-flavour QCD in two dimensions that have 
the quantum numbers of baryons and mesons, constructed out of quark and anti- 
quark solitons. In isolation the latter solitons have infinite energy, corresponding 
to the presence of a string carrying the non-singlet colour flux off to spatial infinity. 
When N c solitons of this type are combined, a static, finite-energy, colour singlet 
solution is formed, corresponding to a baryon. Similarly, static meson solutions 
are formed out of a soliton and an anti-soliton of different flavours. The stability 
of the mesons against annihilation is ensured by flavour conservation. The static 
solutions exist only when the fundamental fields of the bosonized Lagrangian belong 
to U(N c xNf) rather than to SU(N C ) x U(Nt). Discussion of flavour symmetry 
breaking requires a careful treatment of the normal ordering ambiguity. Our results 
can be viewed as a derivation of the constituent quark model in QCD2, allowing a 
detailed study of constituent mass generation and of the heavy quark symmetry. 
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1. Introduction 



One of the key outstanding problems in strong interaction physics is the deriva- 
tion of hadron spectroscopy from QCD, the underlying theory. Quarks were first 
postulated as constituents of hadrons to describe qualitatively the spectroscopy of 
mesons and baryons containing the three lightest u, d and s quark flavours! 11 Sub- 
sequently, it was realized that the short-distance properties of strongly-interacting 
matter could be described exactly in terms of current quarks and the asymptotic 
freedom of QCD. The phenomenological successes of current algebra and chiral 
symmetry 131 implied that the current light (u, d, s) quarks must be much lighter 
than the original constituent quarks, and the relation between current and con- 
stituent light quarks awaits clarification.' 4 81 The distinction between current and 
constituent heavy quarks (c, b, t) is also not clear, but is not so crucial. QCD- 
inspired potential models work well for mesons made out of heavy quarks, and lat- 
tice techniques provide fair understanding of the effective heavy quark-anti-quark 
potential. 

The mystery of the relationship between light current and apparently the heav- 
ier constituent quarks is only deepened by the successes of calculations of baryon 
properties made using the Skyrme model 19 ' 101 , a soliton in the low-energy chiral ap- 
proximation' 111 to QCD in terms of bosonic matrix variables. Constituent quarks 
do not appear in the Skyrme model, their roles being usurped by coherent states 
of current quarks. 

In 1 + 1 dimensions, the spectrum and interactions of mesons in QCD2 were 
first discussed in the framework of the large- N c expansion.' 12 1j1 For baryons, non- 
Abelian bosonization methods' 161 applied to QCD2' 171 have made it possible to ob- 
tain the low-lying spectrum' 18 211 in the case of an unbroken light flavour symmetry, 
again without any reference to the idea of constituent quarks* More recently, ex- 
plicit asymptotic static soliton solutions of the bosonized heavy quark theory have 

★ The one-flavour case in the discretized light-cone formalism is discussed in ref. [22]. 
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been exhibited. These have the quantum numbers of quarks and an infinite energy 
associated with a colour flux tube of infinite length.^ 

In this paper we extend this approach by exhibiting static soliton solutions of 
QCD2 that have the quantum numbers of baryons and mesons. These new solutions 
are colour singlets and have finite energy. The solutions with baryon number zero 
are bound states of the quark and anti-quark solitons, while those with non-zero 
baryon number are bound states of N c quark solitons, corresponding to mesons and 
baryons, respectively. They provide a theoretical laboratory in which the concept 
of a constituent quark can be dissected. They also provide insight into the QCD 
description of heavy-light Qq mesons such as the D and B, and baryons with one 
or two heavy quarks, to which the previous heavy-quark effective potential and 
light-quark chiral approaches have not been applicable. We show that the D and 
B mesons are likely to contain OZI-evading densities of quark-ant i- quark pairs that 
are absent in the naive constituent quark description, and could play observable 
roles in their dynamics and decays. 

The paper is organized as follows. In section 2 we show how mesons emerge 
as electron-positron solitons in QED2, developing intuition for the QCD2 case 
discussed in section 3. Section 4 discusses explicit meson and baryon solutions 
in QCD2, in terms of quark and anti-quark solitons. Such solutions turn out 
to exist only when the fundamental fields of the bosonized Lagrangian belong to 
U (N c xNf) rather than to SU (N c ) x U (Nf). Section 5 discusses their semi-classical 
quantization. Section 6 contains comments on D- and 5-meson physics, and section 
7 is a summary and outlook. Formal aspects of mass splitting and normal ordering 
ambiguities are discussed in an Appendix. 



f There are also qualitative and group-theoretical indications that such a mechanism could 
be responsible for appearance of constituent quarks in QCD4, but the relevant dynamics is 
as yet unknown. 151 
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2. Mesons from solitons in QED2 



Some of the interesting nonperturbative phenomena in QCD2 have close ana- 
logues in QED2 and are easy to derive, once the bosonized form of the Lagrangian 
is known. In this section we present a rederivation of the relevant results obtained 
long ago, via Abelian bosonization, by Coleman , and add some new results of 
our own, namely explicit solutions in the case of broken flavour symmetry. We 
believe the reader will find this section useful for developing physical intuition for 
the discussion of QCD2 to follow in the next section. 

The Lagrangian of multi-flavour massive QED in two dimensions is 



(1) 



where k is the flavour index, is the gauge potential, and 

= d^A v - d v An = E^F (2) 

The bosonized version of (1) reads' 231 

1 X c C 

C=- Y^^xk? + t^F 2 + A=F Y, Xk + J2 -mkHkN^ cos V^ X k 



Al^k = — )=^d v Xk 

c (3) 
i>ki>k = — VkN^k cos v^x/t 

71 

OO 

q= dx y^tpi^k = -7^y2[xk(°o) -xfc(-oo)] 
J k ^k 

where Q is the total electric charge, C = 1/2 e 7 ~ 0.891 (7 is Euler's constant), 
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N^ k stands for normal ordering with respect to a scale and we take the 9 

[23] 

vacuum angle to be zero, for simplicity. 

We integrate out F and take = m^, and then "absorb" a factor of C/tc in 
a redefinition of m^. Then the Lagrangian (3) becomes 

C = \ E^) 2 - f^QC Xk? + J2 m l cos (4) 

k k k 

The equations of motion in the static case read 

Xk - Xl "> ~ sin v^Xfc = (5) 

l 

e 2 

where a = — . The static potential is given by 
An 

V = 2a(J2xif + J2 m2 k^ ~ cos ^Xk) (6) 
I k 

Multiplying (5) by x'ki summing over k and integrating, we get 

k 

We are looking for finite energy static solutions, therefore the gradient terms must 
vanish asymptotically, 

X*(±oo) = (8) 

so that 

V(oo) = V(-oo) (9) 

* For a detailed discussion of normal ordering in the presence of different masses, see the 
Appendix. 
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Taking Xk(— 00 ) — 0, we see that V(— 00) = 0, and hence V(oo) = 0. Thus 

1 (10) 
cos V^47n\^(oo) = 1 

We thus see that only states with zero total charge Q are allowed, which is what 
one expects, since QED2 is confining. From (10) it follows 

Xk(°°) = Vnn k , n k = 0,±1,±2,... (11) 

For two flavours, from (10) and (11), taking the lowest non-trivial n^'s, namely 
n\ — l,n2 — —1, we obtain 

Xi(oo) = v 7 ^ 

(12) 

X2(oo) = -y/n 

The boundary conditions (12) correspond to a meson built out of a soliton and an 
anti-soliton. Equations (5) with boundary conditions (12) can be solved explicitly. 
When mi = 1712 it is easy to see that xi( x ) = ~X2( X ) and the "string tension" term 
proportional to a in (5) vanishes, leading to two "mirror" decoupled sine-Gordon 
equations for xi, X2- When mi 7^ TO2, a solution can be found numerically. It is 
particularly interesting to examine the solutions for widely unequal masses. Some 
representative solutions are shown in Fig. 1. It is worthwhile pointing out that for 
a mf, m\ the widths of xi and X2 are governed by \ jm\ and l/m2, respectively. 
The lighter "quark", say %2 ? will be more influenced by the string tension, getting 
heavier with increasing ajm\. This is just the intuitive picture we usually associate 
with a light current quark getting a constituent mass of the order of the gauge scale. 

Let us now consider the case m\ — > 00. Then, in (5) with k — 1, we can neglect 
the a term and get a free soliton of mass m\, which in the limiting case tends to 
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a theta function 

Xi 0*0 — (13) 

mi— >oo 

as follows from the explicit form of the sine-Gordon solution, 

Xi = —j= tan -1 exp (\/~Anm\x\ (14) 
V 71 " ^ ' 

Then, for the light flavour, k = 2, 

X2 - 4«X2 - V^Trml sin v^Xi = 4a v / 7r6 l (a;) . (15) 

We see that the light "anti-quark" field X2 feels a point-like "source" term due to 
the heavy "quark" xi- m the following sections we shall demonstrate that a very- 
similar phenomenon occurs in QCD2. 

3. Hadronic Solitons 

Two non-Abelian bosonizations of QCD2 have been developed, one in terms of 

[18] 

SU(N C ) x U(Nf) bosonic variables where N c is the number of colours and Nf is 
the number of flavours, and the other' 191 in terms of U(N c xNf) bosonic variables* 
We will now show that in the U(N c xNf) scheme, but not in the SU(N C ) x U(Nf) 
scheme, there are static solutions that have the quantum numbers of baryons and 
mesons, constructed out of quark and anti-quark solitons. In isolation the latter 
solitons have infinite energy, corresponding to the presence of a string carrying the 
non-singlet colour flux off to spatial infinity. When N c solitons of this type are 
combined, a static, finite-energy, colour singlet solution is formed, corresponding 
to a baryon. Similarly, static meson solutions are formed out of a soliton and an 
anti-soliton of different flavours. The stability of the mesons against annihilation 
is ensured by flavour conservation. 



* The specific case of SU(N C ), Nf = 2 has also been considered in a mixed Abelian non- 
Abelian formalism in Rcf. 24. 
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The effective exact U(N c xNf) ffavour symmetric bosonic action with gauge 



fields integrated out reads 

o2 



[19] 



S eff [u] = S [u] + d 2 xTr dZ 1 (ud-u^ +m' 2 N jh d 2 xTr (u + u^ 



(16) 



where 



[16] 



S [ u ] = l- [ d 2 xTi(dr l u)(d fl u^ + -l- [ d 3 xTrj jk (uidiu)(uidju)(u%u) (17) 
8n J • 12n J 

B 

e c is the strong coupling constant (having dimensions of mass in 1+1 space time), 
u is an element in U(N c xNf), the subscript c denotes projection onto the colour 
part, i.e. averaging over flavour and subtracting the U(l) part, and 

_ C 



/2 

m = m q m — 

Z7T 



(18) 



where m is a normal ordering scale, and Nfh stands for normal ordering with respect 
to in (see the discussion in the Appendix). We shall look for solutions with u in 
diagonal form (see the discussion in ref. [18]), 



^aa'jj' — $aa' ^jj'G 



%y 4-7TX 



(19) 



a = l,...,N c 
j = l,...,N f 

The Wess-Zumino term vanishes for either static or diagonal solutions. In general, 
for non-equal masses, in the static case, 

-i 2 



aj J a J 

aj 



Xal Xf3l 



(20) 



where a c = e 2 /A-K and rrij stands for the j-th mass. The relation of mfs to the 
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"current" quark masses and to a c , as obtained through the normal ordering, is 
discussed in the Appendix. 



The Hamiltonian density is 



a.) 



V = 2a c J2 



^^Xal X/31 



(21) 



+ l] m j 1- COS V^Xaj] 



aj 



where we have added a constant to V, to make V vanish for x a j = 0- The 
equations of motion are 



Xaj ~ 4 "c 



Xal X(3l 



- Vl^m] sin V4^x a j = (22) 



To obtain an integral of motion, we multiply by x' aj an d sum over a,j, to obtain 
in the static case 



E-U x ^) 2+2ac S 



aj 



-i 2 



^^Xa£ X.(3l 



f3l 



+^^m| 1 — cos V^nXaj = (const. 



aj 



(23) 

Namely, the left-hand side of (23) is independent of x. For static solutions, the 
existence of such an ^-independent integral of motion is the analogue of energy 
conservation for solutions evolving in time. 

If we choose the boundary condition Xaj(~ oo) = 0, then the constant on the 
right-hand side of (23) vanishes, and hence also for Xaj (+°o) (denoted hereafter 
simply by Xaj), we must have 

-i 2 



(31 



+ ^ mj 1 — cos a/47tx 



ay 



(24) 



+00 



Note that this condition is also obtained by requiring finite energy, as follows from 
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eq. (21). We infer from (24) that 



1 

^Xaj = n a j integers 

(25) 

^Xai — V^^ n al — V^n independent of a 
l l 

The baryon number of any given flavour I is given by 

Bi = J2 n *i ( 26 ) 

a 

Combining eqs. (25) and (26) we get the total baryon number 

B = Y,Bi = nN c (27) 
l 

which clearly is an integer multiple of N c . 



4. Explicit solutions and bosonization schemes 

For the mesonic solutions, we need B = and hence n — 0. Let us consider 
the case of a meson containing a quark of flavour I = 1 and an anti-quark of flavour 
I = 2, and no other constituents. Thus, we need 

^2 n al = 

I 

(28) 




Note that in the product scheme, i.e. with h in SU(N C ) and g in U(Nf), we would 
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have 

rial =Pa + Ql 
V n (29) 

Then 

Nf Pa + J2<ll = ( 3 °) 
l 

which implies that p a is independent of a and hence must be zero. Thus 

n a i = qi integers (31) 

l 

( I 1 = 1 (32) 
N c qi= l-l 1 = 2 
{ I > 3 

which is impossible. Thus there are no solutions with mesons as quark-anti-quark 
states in the "product scheme" of SU(N C ) x U(Nf). 

Let us now give an example of a solution in the U(N c xNf) scheme. In view of 
the discussion above, such a solution will obviously have components in what we 
called / in ref. [19], namely a non-factorizable part in the colour-flavour space. 

As obtained above, the asymptotic boundary conditions are: 

Xaj(-oo) = 



Xaj ( + 00) = v 7 ™. 



(33) 



aj 



where the set {n a j} must satisfy the constraints (28). A possible solution is 

nn = 1 

(34) 

7112 = "I 

with all other n a j being zero. Having specified the asymptotic boundary conditions 
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at x — > ±00, we must now see whether a solution exists for all x. In general such 
solutions can only be found numerically* The case of an exact SU(Nf) symmetry, 
i.e. equal quark masses, is an exception where an explicit analytical solution can 
be found: 

2 



XllW = - XI2W = ^tan-'[exp(^ m ,)], (35) 



with all others vanishing identically. In order to show that this is a solution, we note 
that with (35) the string tension term proportional to a c in (22) vanishes identically, 
and the equations for the various Xaj decouple, yielding a set of independent sine- 
Gordon equations with equal masses. It is then obvious that with the boundary 
conditions (33), (34), xu( x ) * s the usual sine-Gordon soliton. 

Recently there has been much discussion in the literature of the so-called heavy 
quark symmetry 126 311 . Here this symmetry manifests itself in a rather clear way. 
Consider a Qq meson made out of a heavy quark Q and a light anti-quark q, such 
as the D- or B-mesons. When ijiq is much larger than the scale of the theory, 
m Q ^> e c, its profile tends to a theta function, 



XiQ — ► V^O(x) (36) 

m Ql e c — * 00 



as in the QED2 case, while the baryonic current of the heavy quark tends to a 
delta function, 

J% = A=Y,d x XaQ — ► (37) 
V 71 " — mq/e c ^oo 



Q 

a 



Thus the heavy quark acts as a static colour source, while the profile Xaq of the 
light quark becomes independent of ijiq, as shown in Fig. 2. The physics of B and 
D mesons in the context of the heavy quark symmetry will be discussed in more 
detail in section 6. 



[25] 

★ We have used the subroutine package COLSYS. 
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The presence of the static colour source (37) makes the total energy of the 
Qq system finite. From (21), (26), (33) and (34) we see that had there been no 
other colour source, the light flavour profile xiq on its own would correspond to 
a configuration with a net baryon number —1, one "unit" of colour charge and a 
finite energy density per unit length, associated with a colour flux tube of infinite 
length, 

Xlq — ► 

X — > oo 

x N 2 (38) 



V- q {x) — 2na c [l- w 
x—>oo y JM C 

resulting in the total energy being infinite. This is reminiscent of the single quark 
solution discussed in ref. [8]. It gives precise meaning to the intuitive concept of 
quark confinement: isolated quarks have infinite energy because flux conservation 
forces them to emit a flux tube which has no sink to absorb it. 

Multi-quark baryonic solutions can be obtained in a similar way to the meson 
solutions. For example, taking N c = 3 and B = 3 (a 3-quark state*) we find n — 1 
(cf. eq. (27)). One possible solution is 

n\i = ri2\ = n 3 i = 1 (39) 

Corresponding to a "•uW-like baryon, i.e. a "A ++ ". When quark masses are 
equal, mi = mi = ... = m,N f , this solution coincides with one found earlier 

[18] 

in the "product scheme" , in which the colour part is "frozen", i.e. the string 
tension vanishes identically. The vanishing of string tension in (39) is caused by 
a mechanism similar to the one occurring in the meson with equal quark masses 
described above, leading to 

XliO^) = X2l(x) = X3l(x) = —j= tan -1 exp (v^mrr) , (40) 

This baryon solution is plotted in Fig. 3(a). 



* In our normalization a single quark carries one unit of baryon number. 
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There is an important difference, however, between the meson and the baryon 
cases. The vanishing of the string tension term in the meson is a manifestation 
of the fact that the chromoelectric fluxes of the quark and anti-quark cancel each 

[23] 

other. This phenomenon has its counterpart in QED2 , as discussed in section 2. 
The cancellation of fluxes of N c quarks has no such counterpart and can only occur 
in a non-Abelian theory. Another difference between the baryonic and mesonic 
solitons is that the latter do not exist in the "product scheme" . 

When the boundary conditions (39) are taken together with non-equal quark 
masses, we still obtain a u uuu" -like baryon, with the classical solution (40), as 
in the "product scheme". The quantum fluctuations, however, must be treated 
differently, due to flavour symmetry breaking (see the Appendix for additional 
discussion). 

An intrinsically new 3-quark solution occurs when more than one flavour ap- 
pears in the nontrivial solution, for example, 

nn = n 2 2 = n 33 = 1 (41) 

This corresponds to a baryon in which each quark has a different flavour. Such a 
solution is particularly interesting when quark masses are taken to be non-equal, 
corresponding to a "uds"-like baryon, or to a baryon in which one quark is much 
heavier than the QCD scale, such as the S c or again serving as a theoretical 
laboratory for the study of the heavy quark symmetry discussed earlier?" 

One can also study baryons containing two heavy quarks, as shown in Fig. 3(b). 
For N c = 3, the light quark distribution in such a baryon is the same as in a Qq 
meson. The physical reason for this is that the two heavy quarks are essentially 
at rest and act as a static colour source. They are colour triplets and combine to 

f In 1+1 dimensions there is no spin and therefore baryons must be in the completely sym- 
metric representation of SU(Nf), as the space part of the wave function is symmetric for 
lowest energy states. Thus with N c — 3 and three light flavours there is only decuplet and 
no octet, and there is no analogue of A c . 

I For a different approach to heavy-quark baryons in the chiral soliton framework, see ref. 32. 
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a colour anti-triplet, 3 <8> 3 — > 3*, i.e. the effective field seen by the light quark q 
in a QQq baryon is that of a static anti-quark. This gives precise meaning to the 
concept of constituent quark in QCD2. The QQq case is particularly clear, since 
this type of baryon contains only one light constituent quark, while in Qqq or qqq 
baryons there are two or three such objects, superimposed nonlinearly. A solution 
of the QQq type with mi <C mi <C mj, is shown in Fig. 3(c). 

5. Semi-classical quantization 

As has already been mentioned, one of the specific applications which moti- 
vated this study was that to mesons containing just one heavy quark Q. In the 
constituent quark language, these would be described as Qq mesons, where q rep- 
resents some light constituent quark (u, d or s). Examples include the D, D s and B 
mesons. To describe such mesons within our QCD2 approach, we need to consider 
a quark mass matrix with one heavy eigenvalue M ^> e c , and Nf — 1 (typically 
three) light eigenvalues m <C e c . In such a case, the light quark degrees of freedom 
should be quantized semi-classically, as was already done for baryons made out 

[18] 

of light current quarks. The resulting lump is the best QCD model we can de- 
rive for the concept of a light constituent quark. However, clearly it is a coherent 
state containing in some sense an infinite number of current quarks, at least in the 
massless limit. 

When discussing quantum fluctuations around the classical solution u c (x), we 
write 

U = Au c {x)A- 1 (42) 
with A unitary. Now, we take A to be a product 

A = A c Af (43) 

where A c acts only in the colour space and Af in flavour space. This factorization is 
implied by the fact that although the bosonic variables u(x) belong to U(N c xNf), 
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the full theory is not U(N c xNf) invariant. Only the free-fermion theory has the 
full U(N c xNf) symmetry, and gauging of the colour part results in this symmetry 
being broken down to SU (N c ) x U (Nf) x U(l). If we take A c to be a function of x_ 
only, we see that it will not change the interaction term in S e jj[u] of eq. (16), thus 
having no a c effects. Thus we use "light-cone quantization" for the colour part. 
For Af we consider a function of t only, thus reducing the treatment of quantum 
fluctuations to that which has already been treated in ref. [18], for the equal mass 
case. The unequal mass case is discussed in the Appendix. 

Taking the colour part A c to be a function of X- , while the flavour part A f 
is taken to be a function of t, might turn out to be problematic. A more natural 
choice would have been to take both A c and Af as functions of t, since we are 
performing the quantization around a static classical solution, u c = u c (x). 

In any case, we expect that for equal quark masses the fluctuations in the 
colour space will not influence ratios such as flavour content of hadrons (sec. 6), 
since the physical states are colour singlets. The situation might change, though, 
when masses are different. The flavour content will then in general depend on 
the mass ratios, and the different masses could in general be influenced differently 
by the colour fluctuations. This is because the flavours which are light compared 
with the typical scale of the theory are in the strong-coupling regime, and will be 
more affected than the flavours which are much heavier than the dynamical scale 
of the theory. In our case, we will be interested in the u, d and s contents in the 
approximation where they are all very light. 
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6. Comments on D and B physics 



In the previous sections we have described QCD2 solitons which could serve as 
models for Qq mesons such as the D or B. In addition to giving some insight into 
the concept of a constituent quark from the point of QCD, this study may also give 
some new insights into the dynamics and weak decays of D and B mesons. In par- 
ticular, we would like to comment on the existence and possible phenomenological 
role of non-valence quarks in the D and B meson wave functions. 

There are various phenomenological indications that the proton wave function 
contains a significant density of non-valence ss quarks. Moreover, their abundance 
relative to uu and dd quarks is qualitatively reproduced by Skyrme model soliton 

calculations. The relative abundances of ss, uu and dd have also been calculated 

[21] 

in baryonic QCD2 solitons, and found to be qualitatively similar to the QCD4 

[33—35] 

results. In this chapter we invert the logic: QCD2 mesonic solitons contain 
calculable non-valence quark densities, and we would expect QCD4 mesonic solitons 
and hence physical D and B mesons to contain similar non-valence quark densities. 

The calculation of the different light quark densities in a QCD2 mesonic soliton 
parallels very closely that in a baryonic soliton: [21] 



where <C 3> denotes an average over the collective coordinate representation of 
the mesonic soliton. In the physical case with one heavy and three light quarks we 
find for a Qq\ meson: 



This parallels the previous result for the sea- and valence-quark content of 
baryons in the semi-classical quantization of QCD2 with iVj light flavours and 

* Qq mesons in QCD 2 have also been recently studied in the large- N c limit! 361 



(M\q iqi \M) = < z*Zj > 
(M\q jqj \M) «^» 



(44) 



(M\q iqi \M) : (M\q 2 q 2 , qmW) = 2:1 



(45) 
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N c colours. For a baryon B containing k light valence quarks of flavour v 




k + 1 



(46) 



N t + N c 



and 



((qq)sea) B 



1 



(47) 



where (qq) se a refers to the non- valence quarks in the baryon B, and the total 
flavour content is normalized to f. Equations (46) and (47) suggest an "equipar- 
tition" for valence and sea, each with a content of 1/{N[ + N c ). Since the semi- 
classical quantization of light flavours in Qq mesons parallels that in the baryons, 
the "equipartition" result applies to the light flavour content of Qq mesons as well. 
This implies, for example, that 



and similarly for B mesons. We might expect the ratio (48) to be qualitatively 

[21] 

similar, though possibly smaller by about a factor 2, for the realistic case of 
QCD 4 with light flavour SU(3) breaking! 33 " 351 

The presence of significant amounts of ss quarks in the D and B mesons implies 

[37] 

that the naive OZI rule forbidding disconnected quark diagrams can be evaded. 
This might also have implications for D and B production dynamics, but here we 
only emphasize some possible implications for D and B decays. 

1. Annihilation diagrams: cs — > ud could be more important for D° and D + 
mesons than is normally supposed when only the D s wave function is assumed 
to contain strange quarks. 

2. The final states from D and B decays could contain more strange particles 
than is normally supposed when ss pairs need to pop out of the QCD vacuum 
or be created at the weak vertex. This could help explain the surprisingly 



(D\ss\D) / (D\uu\D) 



1 



(48) 



2 



17 



large [38 ' 39,401 branching ratios for D° — > (j)K° and D + — > <f)K + . This observa- 
tion could also have implications for attempts to estimate the ratio|V^s/V^| 
of Kobayashi-Maskawa matrix elements on the basis of strange final states.' 411 
It might also have implications for the ratios of D — > KK and irir final states. 

3. The ss pairs could provide an additional source for B — > + X. At this time 
it is premature to compare this with the data. 

Detailed investigations of these possibilities should await more realistic calcu- 
lations in QCD4, however. 

7. Summary and outlook 

We have shown in this paper that the spectrum of QCD2 includes finite-energy 
mesonic and baryonic solitons which contain one heavy quark. These solitons 
can be regarded as bound states of the infinite-energy single-quark solitons that 
we found previously. They provide meaning for the previously fuzzy concept of 
a constituent quark, at least in QCD2. A particularly interesting application is 
to the study of mesons and baryons containing both heavy and light quarks. A 
constituent light quark is seen to be a semi-classical coherent state containing 
an indefinite number of light qq pairs, among which non-valence flavours have as 
much as one half of the abundance of the valence flavour. This observation could 
have phenomenological implications for the dynamics and weak decays of D and 
B mesons. 

The next step in the programme of developing accurate QCD descriptions 
of constituent quarks and Qq mesons is to extend the analysis of this paper to 
the realistic case of four dimensions. This may be possible for the lowest-lying Qq 
mesons if they are describable by spherically symmetric wave functions, which could 
be analyzed using an effective two-dimensional field theory in (r, t) coordinates. We 
are now investigating this possibility. 
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APPENDIX 

We discuss here the more formal aspects of mass splitting. In the following we 
will focus mainly on the bosonization in the product scheme, since the technicalities 
are simpler there, while the conclusions are the same as in the U(N c xNf) scheme. 

MB] 

We recall the initial bosonic action in the product scheme 



S [g, h, Ay] =N C S [g] + N f S [h, M~^J d 2 x Tr V" 
+ CmNrn J d 2 x{Tr (M g) Tr h + h.c] 



(A.l) 



where g G U(Nf), h G SU(N C ), is colour gauge field, M = diag(mi, . . . ,m,N f ) 
is the mass matrix, mi are the quark masses, iV™ stands for normal ordering with 
respect to rh, and is the usual non-Abelian field-strength tensor. Now define 
the "dimensionless mass matrix" D as 



D = —M (A.2) 
mo 



where mo is an arbitrary mass parameter. In the gauge A- = 0, integrating out 
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A+, and then taking the strong-coupling limit, e c all rrii, we get 



[18] 



S [g] = N C S [g] + m 2 N m / d 2 xTr \d{ 9 + g f ) 



(A.3) 



m 



CN c m 



elNf_ 
2tt 



i 



(A.4) 



P 



(A.5) 



N C (N C + N f ) 

In the U(N c xNf) bosonization scheme the mass term, denoted M.T., has the form 



M.T. = moCmNm Tr 



D(u + v)) 



(A.6) 



where u G U(N c xNf) and D is the dimensionless mass matrix (in flavour only). 
In the strong-coupling limit the mass term (A.6) exactly coincides with that of 
(A. 3). [19] The classical low-lying solutions to the action (A.3) are given by the 
time-independent form 



9 = diag 



exp (^ i \jjj-^ x ) S j , • • • ,exp (^i\jjj-<i>N f {x) 



(A.7) 



yielding an action density 



s d [g] = - dxJ2 



Ni 



i=l 



2 V dx 1 




<Pi-l 



(A.8) 



with mass parameters 



Each <pi has a classical solution 



rn 



ttia 



1 ™2 



m 



m 



(A.9) 



<f>i(x) 



4Nr 



arctan 



7T 



8ir 

<>X P 1 \IJ^ m i X 



(A.10) 
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with an energy 



Ej = 4m, 



2N r 



TV 



(A.H) 



The minimum classical energy solution is given by a solution with exactly one 
non-trivial entry 



g (x) = diag 



1,1,..., exp 



/ 4:71 



(A.12) 



with mjVy chosen to be the smallest mass. 

In order to quantize the system we introduce the collective coordinates A(t), 



g(x, t) = A(t) g (x) A*(t), A(t) £ U(N f ) 



(A.13) 



Writing A 

get 



(A\, . . . , AjVy-i, z) where A\, . . . , Ajy f _i, z are column vectors, 



wc 



§ [g] ~ ~ S N = 2~M N ~ f J dt ( Dz )i( D z)i + f dt&l 



2tv 



M Nf N c 



N f 



(A.14) 



2 ~ 2 m |2 

m i -"ijv>)l^i| 



.Dz — z — z{z^ z) 



Nr. 



2Mj 2ir 



l 4av 



y/2 (Nc 



1 — cos \ — 4>i 1 dx — — , 

N c J mi \ TV 



i = l...Nt 



The resulting Hamiltonian is 



H = Am^ [ 



2N r 



TV 



1 + 



TV 



2N, 



C2(R) - §-{N f - 1) 



N f 



mi - m N 2 

\Zi\ 



i=i 



(A.15) 

The Hamiltonian depends on uiq only through m/Vy, therefore the overall mass 
scale is undetermined, and only mass ratios are meaningful. The quantity \z{\ 2 
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is proportional to the qiqi operator. Therefore the last term in (A. 15), due to 
mi 7^ Mjq f , is proportional to the weighted average over flavours of the qq operators, 
the weights being proportional to the mass differences. Since this term comes from 
quantum fluctuations around the classical solution, consistency with the semi- 
classical approximation requires that it be very small compared to one. 

Generalizing eq. (A. 2), by introducing an extra undetermined mass parameter 
rriio for each flavour, we define pi = rrn/miQ. The mass term has the form 

N f 

M.T. = CmNrn Tr (Mg)=Cj2 momN A p igii (A. 16) 

i=i 

By integrating out the gauge fields and taking the strong-coupling limit, we get 

S [g] = N C S [g] + £ m?^ / d 2 x^(g tt + (A.17) 

J mo 



p+i 

(A.18) 

The classical low mass solutions to the action (A.17) are given by (A. 7) with an 
action density (A. 8), but now with mass parameters 

~2 - 2 m i ~ 2 / a i rVi 

m i = Pim i = m, (A. 19) 

mo 

Each (pi has a classical solution (A. 10) with energy (A. 11). The minimum classical 
energy solution is given by the ansatz (A. 12), with rnjV/ chosen such that it is the 
smallest mass. Define now A(t) as in (A. 13), resulting in an effective action like in 



m. 



CN c mio 



2tt 
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(A. 14), and a Hamiltonian 



H = M<1 + 



7T 

2N r 



C2(R) - §- f (N f - 1) 



M = 4m Nf 



N f 

+E 

i=l 



//// / /»/o \ 1+p 



m Nf \m N0 J 



7T 



- 1 

(A.20) 



There are several problems connected with these arbitrary mass parameters. Since 
the choice of the parameters m^o is arbitrary, we can choose 



mo = m N0 



mi 



m N . 



p+i 

p-i 



in which case H has no explicit contribution due to the mass differences, which is 
unacceptable. 

The question now arises at what stage should one apply the semi-classical 
approximation.' 201 The classical solution is governed by rhi but the lightest rhi 
is not necessarily the lightest mj. In fact, the normal ordering ambiguity is 
present even in a non- interacting theory. To see this, we formally set N c — 1, 
i.e. we have the field g with level 1 which means free quarks. There is no gauge 
group, so the dependence on e c should vanish and indeed p = 0, rhi = CrriiQ and 
/ 2m Nf m N0 \ 



M = AC 



V 



7T 



1 



. The lowest multiplet corresponds to a Young tableaux of 



one box and C2(R) = 



N 



f 



1 



2N f 



-, therefore 



H = M 



7T\ 2 Nf 

2, 



N f 



i=i 



m Nf m N0 



(A.21) 



which in general does not look like a non-interacting theory, unless one takes m^o = 
mj\fo for all i. 
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One could argue that the above example is not relevant for a generic case, 
since for small N c the semi-classical approximation is not justified. This can be 
seen from (A. 15), as follows. Consistency with the semi-classical approximation 
requires the quantum corrections to the energy to be much smaller than the classical 

contribution. In the case of one box this means (S ) — <C 1, which for any 

Nf > 2 is not correct. 

Nevertheless, since the choice of the m^o's is a priori arbitrary, we see that 
there is an ambiguity. In the non-interacting case it is obvious what is the right 
choice of the m^o's, but in general we do not know what is the choice leading to 
the optimal semi-classical approximation, as is best illustrated by the case of equal 
quark masses, which can be made to appear unequal by a suitable choice of the 

TOGO'S. 

One can summarize the situation as follows. When all masses are equal, the 
treatment of ref. [18], in the strong coupling limit, yields an effective Lagrangian 
in terms of flavour degrees of freedom only, with a mass scale that involves the 
coupling a c , as in eq. (A. 4). When the current mass M of some quark is heavy, 
M 2 ^> a c , we expect the constituent mass of that quark to be heavy and close 
to the current mass. The problems arise in the intermediate cases, when we do 
not know what is the "best" starting point for the normal ordering scale, before 
going to the semi-classical approximation. Of course, if one were able to sum all 
corrections, one would obtain the full result, regardless of the starting point. 
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FIGURE CAPTIONS 



1) Mesons in QED2 with two flavours, solutions of eqs. (5), (12). 

(a) a = 0.1, mi = 0.5, m 2 = 0.5; 

(b) a = 0.1, mi = 2.0, m 2 = 0.5; 

(c) a = 0.1, mi = 4.0, m 2 = 0.5; 

(d) a = 2.0, mi = 4.0, m 2 = 0.5. 

In all four plots the upper line corresponds to xi> and the lower line to X2- 
(a), (b) and (c) show how, with increasing mi, the xi profile converges to a 
step function, while the light flavour profile, X2, remains almost unaffected. 
The a dependence can be inferred from the difference between (c) and (d): 
X2 gets heavier with increasing a/m 2 . The effect is not a sharp one, because 
part of the a dependence is already included in the value of m 2 , through the 
normal ordering (see Appendix). 

2) Mesons in QCD 2 with N c = 2, N f = 2, solutions of eqs. (22), (34). 

(a) a c = 1.0, mi = 0.1, m 2 = 0.1; 

(b) a c = 1.0, mi = 0.1, m 2 = 0.5; 

(c) a c = 1.0, mi = 0.1, m 2 = 2.0; 

(d) a c = 10.0, mi = 0.1, m 2 = 2.0. 

X11 and X12 are the upper and lower continuous curves, respectively. The 
"non- valence" components X21 and X22 are denoted by dot-dashed and dashed 
lines, respectively, except for (a), where they are exactly zero. 

3) Baryons in QCD 2 with N c =3, Nf—3, solutions of eq. (22). 

(a) a "wW-like baryon, eq. (39), with a c =1.0, mi=m 2 =m3=0.1; 

(b) a "-uc6"-like baryon, eq. (41), with a c =1.0, mi=0.1, m 2 =0.5, m3=0.8; 

(c) a u ubt" -like baryon, eq. (41), with a c =1.0, mi=0.1, m 2 =1.0, m3=2. 

In (a) the "non- valence" components, X21, X3b etc are exactly zero. In 
(b) and (c) the continuous lines denote the "valence" components, while 
the dashed and dot-dashed lines denote the "non-valence" components. The 
"non-valence" components of the heavy flavours are clearly much more sup- 
pressed than those of the light flavour. 
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